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Abstract 

In order to model the phase-coherent scattering of electrons in two-dimensional electron gases in the 
presence of Rashba spin-orbit coupling, a general partial-wave expansion is developed for scattering from 
a cylindrically symmetric potential. The theory is applied to possible electron flow imaging experiments 
using a moveable scanning probe microscope tip. In such experiments, it is demonstrated theoretically that 
the Rashba spin-orbit coupling can give rise to spin interference effects, even for unpolarized electrons at 
nonzero temperature and no magnetic field. 
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I. INTRODUCTION 



There has been recent interest^*^ in utilizing the spin degree of freedom in semiconductor 
devices, where the charge carrier's spin provides an additional degree of control and flexibility 
towards developing devices that are faster and more efficient devices than conventional electronic 
devices. One component of potential "spintronic" devices, the spin transistor proposed by Datta 
and Das'^, modulates the current passing through a semiconductor due to the presence of the spin- 
orbit interaction, which couples the electron's spin with its kinematical motion. Interest in the 
spin transistor has generated numerous theoretical and experimental investigations into the spin 
dynamics under the spin-orbit interaction in two-dimensional electron gases (2DEG). 

In layered semiconductors devices, the two predominant sources of spin-orbit coupling arise 
from either structure inversion asymmetry (SIA or Rashba interaction^) or bulk inversion asym- 
metry (BIA or Dresselhaus interaction^). The BIA spin-orbit interaction arises from the breaking 
of inversion symmetry by the inherent asymmetry of the atomic arrangement in the structure and 
is not very amenable to external manipulation. The Rashba spin-orbit coupling, on the other hand, 
arises from band bending at the interfaces between semiconductor layers and/or any external elec- 
tric fields applied to the the device. Unlike the Dresselhaus coupling, the strength of the Rashba 
coupling can be partially controlled by application of an external electric field^ and in principle can 
be made the dominant form of spin-orbit interaction in the 2DEG. Such tunability of the Rashba 
interaction is ideally suited for applications in spintronic devices, and as such, only the Rashba 
spin-orbit coupling will be considered in this study. 

Numerous studies have been conducted on the diffusive transport of spins in the presence 
of spin-orbit couplin g^i^i'" in order to investigate a variety of phenomena, such as the spin Hall 
effect' Most of the studies were conducted up to the first-Born approximation for the scat- 
tering from nonmagnetic impurities, and the results were disorder averaged. However, there are 
many cases where such statistical theories are not warranted. For example, coherent scattering 
from a fixed set of impurities, which give rise to quantum interference effects induced by multiple- 
scattering events from the localized impurities, can't be described by such statistical theories. One 
method of tackling such problems is multiple scattering theory, which has been routinely used in 
optical and acoustic scattering and has been proposed as a method for understanding the fringing 
patterns in recent imaging experiments on electron flow in 2DEGi^. In scattering theory, the 
effect of a scatterer k can be localized to a point in space at the center of the scatterer, f^, such that 
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an operator, T^, can be constructed which generates the scattered wave, '^s{R)i from the incident 
wave, ^ia{R), evaluated at the site of the scatterer, R — fk. 

The subscript, f = fk means to operate Tk{R, r) upon $in(r) and evaluate the result at r = rfe. In 
the presence of N point scatterers, the total wave function is then given by 

TV 

^{R) = ^-^{R) + V (fk{R, r)^(r-)) (2) 

Thus the complete wave function can be found if (Tk(R, r)^(r) j is known at each scatter k. 

In the following article, a multiple-scattering theory in the presence of Rashba spin-orbit cou- 
pling in a 2DEG is developed. The general formalism is presented, along with the explicit calcu- 
lation of the scattering operator, Tk, for a cylindrically symmetric well/barrier, which will be used 
as a model for impurities in a 2DEG. As an application, the methodology is applied to possible 
flux measurements for phase-coherent transport in a 2DEG with Rashba spin-orbit interaction in 
the presence of a scanning probe microscope (SPM) tip in zero magnetic field. Additional interfer- 
ence effects arise in the flux measurements due to spin interference effects caused by the Rashba 
coupling. 



II. SCATTERING FROM A CYLINDRICALLY SYMMETRIC POTENTIAL: PARTIAL-WAVE 
EXPANSION 

The Hamiltonian for a 2DEG in the presence of the Rashba spin-orbit interaction and impurities 
is given by 

^ Ho + V{x,y) (3) 

where aj are the Pauli spin matrices, and a is the Rashba spin-orbit coupling constant. The eigen- 
states and corresponding eigenvalues for the free-particle Hamiltonian with Rashba spin-orbit cou- 
pling. Ho, are given by 

|A;(^),±(^)) = Iky = kcos{e),kx = ksm{e))\±{e)) (4) 
E. = '^^ak (5) 



where 



tan(^) = ^ 

Ky 



and k = \Jk\^- ky. 

The dispersion relation in Eq. ^ represents two parabolic bands centered upon k = ±m*a/h^. 
For states propagating with their momentum vectors making an angle 9 with respect to the F-axis 
and for an energy E > 0, there exists a two-fold degeneracy with the degenerate states given by 

IhWMO)} = (7) 

where ki(2) = ki(^2){cos{6)Y + sin(6')X) with 



m*a I m*a\ 2m* E 



The states \ki{6), + {6)) and \k2{6), —{9)) represent plane- wave states whose spin states are quan- 
tized in the plane, perpendicular to the momentum direction. 

In polar coordinates, which are useful when considering scattering from a localized, cylindri- 
cally symmetric potential, Hq can be written as 

i/o = -^f— + -- + -— exp(i^) + 

2m* \dr^ rdr dOy \^ exp(-ie) - ) 

(10) 



The eigenstates of Hq which represent states propagating outward from or towards a particular 
origin, fi, can be written as 

- ^ ^ - a/A?7/ H;^(ki\R-ri\) 

2^2 \-iH{l^{ki\R- ri\) exp{-i' 

- ^ ^ - v^/ Htik2\R-ri\) 



ri|jexp(^-it 



(11) 



FIG. 1: Scattering of an incident plane wave ^^{R) [Eq. ( [751 1 from a potential located at f/., Vk{R) 
[Eq. (fT4b l. The various angles and vectors used in Eq. (flSl and Eq. (fT6l are illustrated. 

where H^{z) are Hankel functions given by Hj^{z) = Ji{z) ± iy/(z), and fci and ^2 are given 
in Eq. A similar solution to Eq. (flOb for a cylindrical well has been given befor o^^i^^ . The 
states xfm)iRi E) satisfy a flux orthogonality condition through a circular surface surrounding 
the origin, fi, which is given by 
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2 / ^Octa.E\\\m))m)\J'rAm))W)\)\xU,E)-m)^(^ = —Jl,m6a,, (12) 



where the current operator, J, is given by 

J = J^X + JyY 

— + -T^Y X + — - -ax Y (13) 

m* n J \m* n J 

The states in Eq. dTTT) can be used to generate the scattering operator, [Eq. ([T])], for scattering 
in the presence of Rashba spin-orbit coupling. The following treatment follows closely a previous 
treatment for constructing in the absence of spin-orbit coupling.^' We will begin by solving the 
Schrodinger equation for an eigenstate of Hq [Eq. ©] incident upon a cylindrically symmetric 
potential centered at r^, Vk{r), as shown in Figure [T] Vk{r) is given by 

Vn for |r — fj,\<a 
Vkir-) = \ (14) 

for |r — r^l > a 
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where a is the radius of the scattering potential. Consider the case of an incident wave propagating 
with momentum k± = k±{cos{6o)Y + sin(6'o)X), where ± denotes a particular eigenstate of Hq 

[\k+{eo),+{0o)) = \kM,+ido)) and |L(^o), -(^^o)) = |4(^o), -(^o)) given in Eqs. ©-(H]. 
The incident wave function, ^^(R), written in a coordinate system centered about the scatterer at 
rfc is given by 



1 - - / 1 

— exp(ifc± ■ fk) exp(ifc±r^ cos(^o - O^J) (^^ ^^p(^_;0^^ 



(15) 



where rg ^ = \R — rA h the distance measured from the center of the scatterer and 0?- is the 
angle with respect to the F-axis of the vector r, i.e., 



exp(iqj = (16) 



The wave function outside of the scatterer, \R — rk\ > a, can therefore be written as 

^tiR) = $±(^) + vl>|(^) (17) 
where the scattered wave function, ^^(R), can be written as 

oo 

^siR) = E f^'xtM^ E) + f^\UR. E) (18) 

The wave function inside the cylindrical potential, {R) , can be similarly written for \ R—rk\ < a 
as 



l = ~OD 

+ [xtM E-Vo) + xliiR, E - Vo)) (19) 

Note that Eq. (fT9l) contains both incoming and outgoing states [as given in Eq. dTTTl l in order to 
remove the Yj terms in which are singular at R = fk. 

From the continuity equations of the Schrodinger equations, ^f{R) and 'i>^{R) must satisfy 
the following conditions for all R such that \R — f^] = a: 

^^{\R -fk\=a) = ^±(1^ -rk\= a) (20) 



h'' dentin) ^ ^ exp(ia 



2m* dr 'l«-.l=a y2m*a \^exp(-i^g) 



2m* ar 'l^-^-fcHa l2m*a I exp(-i^|) 




In the following discussion, aji — a, i.e., the spin-orbit coupling strength is the same inside and 
outside the well. 

The solutions for the various coefficients df" and are given in Appendix A for a cylindrical 
well/barrier. In the following, we are interested in studying the wave function away from the 
scatterer, so the relevant coefficients are fj^^ and which can be written for convenience as 



fi = 2rexp{ik±-fk)exp{-iWo)- 



/±2 = 2i' exp(ifc± • ffc) exp(-i/^o) 4^ (22) 

The coefficients, fj^^ and f^"^, depend upon the energy and the form of the potential but do not 
depend upon the initial direction of the incident momentum vector, 9q (different potentials will 
generate a different dependence of the coefficients upon E, etc.). Thus the scattered wave 

I — * 

function, ^^(-R) can be written as 

oo 



'±2 



where 



1 / mirn^rJA} eM^Oo)Hi{hr^.jBl 

2 



+ - 



1 / Hiihr.,)A^ exp{ieo)Hiik,r^.)Bf 



(24) 
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where 

= ft' + P 

Bj = ft' - P (25) 

where Hi = in Eq. (l24b (the + sign will be implicitly assumed for the Hankel function for the 
rest of this paper). 

The operators T^^i in Eq. (E^ and Eq. appear to generate the 1*^ partial wave from the 
incident wave function, evaluated at the site of the scatterer. The only problem with this 
interpretation are the various factors of exp(iZ6'o) occurring in Eqs. (I23t - (l24ll . Since different 
incident waves will possess or be a superposition of different incident momentum directions (i.e., 

in Figure[l]), an additional operator needs to be constructed which generates the various factors 
of exp(i/^^o) from the incident wave with energy E. The operator, Du can be constructed such that 
for any given state of the form |$in) = ci|A;i(6'o), +(6*0)) + C2|/c2(^o), -(^0)), 

= exp(i/0o)$m(i?) (26) 
The operator, Di which satisfies the above equation is given by 

aiPi hiPi+i 

Di= \ I' ' I (27) 



QPi_i diPi 



where 



P. = Ir^^^-i^^l (28) 

with Po = 1- The construction and full expression for Di is given in the Appendix B. 

The operator which generates the scattered wave from the wave incident upon scatterer k, Tk, 
can finally be written as 

00 

f , = ^ i'exp(-i/^g)GfP, (29) 



l=—oo 



where 

G 



2 I iexp(-i^|)if,+i(fcir^,.J mirn,r,) 



2 \-iexp{-ie)Hi^,{k2r,,) m{k,r.,) 
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where tl^, = A^'/a^ + B'jIc^+i, = A^^^-i + B^^^^du tj) = A^^lfai + B^Jc^+i, and tfi = 
—A'l'^fii^i — B^^'^idi. From the values for the various /f '^^ and /f '^^ calculated in Appendix A for 
a cylindrically symmetric barrier/well, it can be shown that t]}^ = tl^^^ and tl\ = tf'_i. Note that 
the form of is the same for all cylindrically symmetric scatterers; the values of the scattering 
amplitudes, f'^-' , depend upon the actual potential used for the cylindrically symmetric scatterer. 



Ill, MULTIPLE SCATTERING THEORY 

For isolated scatterers, the overall wave function at R can be written as 

N 

^{R) = $,„(i?) + J2iTk^mR=r, (31) 
k=l 

Equation (IbTT) indicates that if the value of '^{R) and its derivatives [due to the Pi dependence of 
Tk in Eq. (12^ 1 at each scatterer is known, the entire wave function ^(-R) is completely determined. 
In principle, the values of ^(R) and its derivatives at each scatterer can be found using Eq. (BTT ). 
In practice, it is only practical to calculate the first few derivatives of ^{R) at each scatterer. When 
the size of the scatterer (or in general, the scattering length) is much smaller than the wavelengths, 
i.e., kia <^ 1 and <^ 1, the only significant contribution to comes from the / = term in 
Eq. (l29b . This is analogous to the heavily studied "s"-wave scattering models, and in the following 
discussion, only the / = term in Eq. (l29b will be considered. 



IV. LOW ENERGY SCATTERING LIMIT 

In the limit fci(2)a ^ 1, Eq. (IbTT) can be approximated as 

N 

^{R) = $i„(^)+ J]GS(^)vl/(r-;) (32) 



k=l 



where Gq can be written as 



2 \^ _iexp(-i^^^)i7i(A;2r^,.J ^o(A;2r^,,,. 
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It should be noted that Eq. (l32l) and Eq. (l33t are similar to the Lippmann-Schwinger equation for 
a potential, V{r) comprised of delta- functions: Vs{r) = J2k ^kS{r — Vk), which has been used 
before in previous studies of the spin dynamics in the presence of spin-orbit couplin g^i'^i^^ . From 
the Lippmann-Schwinger equation, the wave function in the presence of Vs{r) is given by 

N 

= + ^ VkG+{R, n, E)^{n) (34) 



k=l 



where G+{R, r^, E) is the Green's function in the presence of Rashba spin-orbit coupling, which 
is given by Eq. ( IC13I ) in Appendix C and is similar in form to Gq(R). The form of Gq{R) would 
be identical to G^{R, r^) if 

|£ = |l (35) 

In general, Eq. is not satisfied, although for ka <^ 1 Eq. (l35t is approximately correct. The 
difference G+{R, fk,E), and Gq{R) can be understood as follows: G+{R, fk, E) propagates the 
scattered wave function from the a delta-function potential, whereas Gq{R) propagates the scat- 
tered wave function from the finite-sized potential, Vfc(f^ [Eq. ([T4b l. which, in the delta function 
limit (a ^ 0,Vo ±oo, nVoo^ = ±Vk), doesn't scatter (i.e., the scattering coefficients, tl g and 
t|o' vanish). 

Far away from the scatterers (A;i(2)r^ ^ 1), Eq. (l33l) can be written as 



R,rk 



(36) 
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where 



ka 
tk 

exp(i50fc) 
exp(i0fc) 



ki + k2 f m*a 
ki — ^2 am* 



2m* 



? I 
rfc,o''fc,ol 



'-kfi'^kfi 
rfc,0''A;,0l 



(37) 



and U{9, (j)) is a rotation operator given by 



U {6, 0) = exp ( i-az j exp (i6'ax) exp ( -12^^ 



(38) 



Eq. (l33t contains a dynamical factor which depends upon the distance from scatterer k multiplied 
by a sum of two rotation operators. In the presence of spin-orbit coupling, the low energy limit 
scattered wave functions possess an "s"-wave character and also a "p"-wave character due to the 
6*? dependence in Eq. (l33t . which vanishes in the limit as a ^ 0. 

In the calculations to be performed, 5tk/tk ~ 0.05 and 5(j)k <^ 1. In this case, Eq. (l33l) can be 
approximately written as [for A;i(2) |-R — ^fe| ^ 1]: 



nr 



exp 1 1 



krR,r-, - I + 



(39) 



which now contains a dynamical, distance-dependent factor times a single rotation operator, 
U{karj^^^,6§J, which corresponds to a rotation by an angle kaVj^^^ about the spin-orbit field 
for propagation along the direction -§^^. In this limit, the total wave function in the presence of 
N scatterers is given by 



Trr, 



exp 1 1 



kr 



IT 

- + 



U{karn,,^,e^^)<I^{f,) (40) 
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From Eq. (l40b . knowing the value of the wave function at each scatterer k, "^{fk), completely 
determines the total wave function ^{R). The various values of ^(rfc) can be found by setting 
R = r^in Eq. ( B^ for each scatterer k. This provides a system of 2N linear equations to solve for 
the various ^(rfc). The resulting system of equations can be expressed in matrix form as 

M$ = J (41) 

where $ and are 2N by 1 matrices where $(2A; - 1) = "^^n), ^{2k) = "^^{n), 4>{2k - 1) = 
(p^^k) , (f){2k) = ^''^(rfc), and M is a 2A^ by 2A^ matrix where M(m, m) = lform= ltom = 2A^, 
and for k, j G [1, A^] and k ^ j, 

M(2A;-l,2j-l) = -{Giin))^^^ 
M{2k,2j) = -{Giin))^^^ 
M{2k-l,2j) = -{Giin))^^^ 
M{2k,2j-1) = -{Giin))^^^ 

(42) 

where Gi{fk) is given in Eq. (l33l) . Once M is specified, \l/ can be found by inverting M as follows: 

$ = M-i0 (43) 



V. APPLICATIONS TO FLUX MEASUREMENTS IN TWO-DIMENSIONAL ELECTRON 
GASES IN THE PRESENCE OF RASHBA SPIN-ORBIT COUPLING 

Recent experiment o ^ ^ have imaged electron flow in a 2DEG by monitoring the changes in 
conductance through a single quantum point contact (QPC) as a moveable SPM tip is scanned 
above the surface of a heterostructure. By applying a negative voltage between the SPM tip and 
the 2DEG, a potential that electrons in the 2DEG can scatter from is generated. As the SPM 
tip is scanned over the surface of the heterostructure, the backscattered current into the QPC is 
monitored, which is proportional to the the electron flow in the probed region in the absence 
of the tip. In this manner, an image of the electron flow is produced. On top of the electron 
flow pattern, interference fringes can also be observed in these experiments, which are due to the 
interference between the single scattering trajectory from the SPM tip with the single scattering 
trajectories from random impurities within the sample. Recently, additional interference effects 
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FIG. 2: Electrons injected by the emitter QPC (represented by <I>in) are backscattered due to random im- 
purities in the 2DEG along with a moveable "scatterer" generated by an SPM above the 2DEG surface. 
Measurements of the backscattered flux into the detector QPC as a function of the SPM tip position, Rt, 
can be used to image the electron flow from the emitter. 

for electrons emitted from a single QPC were also reported for an SPM tip in the presence of a 
fixed reflector gate^, where interference fringes result not only from the interference between the 
single scattering trajectories of the reflector gate and the SPM tip but also from the interference 
of the double scattering trajectories involving the reflector gate and the single scattering trajectory 
from the SPM tip. 

The samples used in the above experimental studies were GaAs/AlGaAs heterostructures, 
which have very low spin-orbit coupling— (a = 3 x lO^^'^ eVm), and the experimen- 
tal results were well described by quantum simulations/calculations in the absence of spin- 
orbit couplin g'^i^'^i'V'i^3^24 _ However, other samples can possess considerably larger spin-orbit 
coupling,^^ such as InAs, which can have a = 4 x 10^^^ eVm. The question therefore arises as to 
what effects or signatures of spin-orbit coupling exist in electron flow imaging experiments using 
a moveable SPM tip. 

Figure |2l shows the setup under consideration. A point source is used to inject electrons into 
the 2DEG (analogous to a QPC). The injected electrons are backscattered by random impurities 
present in the sample and by the potential generated by the SPM tip placed above the surface of 
the 2DEG. The backscattered current into the detection QPC is then measured as a function of the 
SPM tip position, Rt. In the setup shown in Figure El the possibility that the detection QPC can 
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be separate from the emitter QPC is allowed. Such experimental geometries have been used in 
magnetic focussing experiments in the past^^. 

The injected electrons from the emitter are taken to be an unpolarized beam comprised of an 
equal mixture of the cylindrical wave-like states $i(-R) and $2(-R)j which are given by 



MR) 



exp 1 1 



TT 



1 exp ( 1 



TT 



Xo,T(^) + Xo,i(^) 



(44) 



Each state represents current of hk/m* being injected into the emitter, as shown in Figure |2l Far 
away from the source, the states can be approximated as 



MR) 



= exp (ik\R\^U{K\R\,9e) (I 



7l\R 

1 



7t\R\ 



exp (ik\R\)U{k^\R\,ee 



(45) 
(46) 



where 9e is the given with respect to the emitter. Note that when a — > 0, the states $i and $2 
correspond to a spin up and spin down electron being injected from the point source. 
The net current injected into the detector is given by the following formula: 



d^ (Jx(e) sm{9) + Jy^e) cos{9)) 



(47) 



where 



h 



m* I J R=e n 

h 



Jy{e) 



-Im 



J R=e n 



(48) 



and e = e{cos{6)Y + sin(6')X). In the actual experiment/simulation, the current injected into the 
detector is measured as a function of tip position several microns away from the detector. As will 
be discussed later, the interference between the incident waves and the scattered waves can be 
neglected in Eq. (l47t due to thermal averaging, so only the scattered wave function, ^5, needs to 
be considered. Phase-coherent transport is assumed in the application of Eq. (l47t . Additionally, 
since the width of the detector is taken to be negligible, only the current operator of^s evaluated at 
the detector needs to taken into account in Eq. (l47t . The injected current into the detector becomes 



2e — Im 

\m* 



^^l{Rd)ax^s{Rd: 



(49) 



14 



Using the form of the wave function in Eq. (l40b evaluated at the site of the detector, the injected 
current can then be written as a function of tip position and energy as 



+ 



t2 

rk,d h 
Aehk 



'^Hrk)U\k^n,d,et) 



TT m 



V 

- tktj 



cos 



+ cos( 



+ 



■^\fk)U\k^rk,d, et)U{k^rj^d, 0^)-^{fj) exp(i[fc(rj- - rk^) + 0j - 0fc]) + h.c. 



IT h 



1 X 



TkA) + 



+ h.c. 



4e a \^ 



X 



^\fk)U\kark,d, et)aYU (kaVj^d, Of)^if^) exp(il(rj- ^ - rk,d) + (pj - (pk) + h.c. 

(50) 

where the energy dependence of /i(-Rt, E) comes in through the energy dependence of both k and 
tfc. Since all experiments are done at nonzero temperatures, thermal averaging of Eq. dSOb becomes 
necessary. Assuming the injected current is a result of a small potential drop, 5V , over the emitter 
QPC and that the electrons on both sides of the emitter QPC can be described as being free 2DEG, 
the energy weighting function is simply related to the fermi-dirac distribution function, f{E), and 
is given by: 



-f\E)dE 



( E-Ep 

&E exp(^ 



keT 



1 + -P(^))' 



(51) 



where Ep is the Fermi energy. All quantities calculated will be thermally averaged using Eq. (BTt . 
The thermally averaged injected current is then given by 



^i{Rt) = - / fi{Rt,E)f'{E)dE (52) 
Since n{Rt) contains interference terms that go like A;exp(ifcr) (where r is some length related 
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to the various distances in the system), the following integral^l will be useful in evaluating Eq. (l52b : 

poo 

— I kf' (E) exp{ikr)dE = —2i\T exp (ikpr) smh~^{\Tr) {1 — \Trcoth{\Tr)) 
Jo 

+ =^ exp (ikpr^ sinh~^(Arr) (k^pr — 2ArCOth(ATr) ) 
kp ^ ^ 

— ^ exp (ikpr) sinh~^(ATr)A^r (coth^ (At?") + sinh^^ i^Tr)) 
kp 

= exp{ikpr) smh.~^{XTr)g{r^T, Ep) (53) 

where Xp = kpTrkBT{2Ep)~^. For large r, interference terms in fJ,{Rt) decay as r exp(— At^). For 
T = 3K, Ep = 16 meV, m* = 0.022mo (where mo is the free electron mass), Ay = 2.35/im~^, 
which allows one to neglect the interference between the incoming wave and scattered wave when 
calculating fi{Rt) many microns away from the QPC. 



A. The Single-scattering limit 

Before considering the case of multiple-scattering (which can only be analytically solved for 
simple cases), it is useful to consider the single- scattering case for an unpolarized beam (i.e., 
averaged over the incident waves $i [Eq. (l44l) 1 and $2 [Eq. (l44b l). In this case, the value of the 
wave function at scatterer k is given by 

^i(rfe) = —^exp{ikr,^k)U{Kre,k,0'^)r] (54) 

for incident wave $1 [Eq. (l45t l and 

^2(6.0 = —^exp{ikr,^k)U{k^r,^k,0'^)r] (55) 

for incident wave $2 [Eq. (l46b l. 

First consider the case when the detector and the emitter are one and the same. The spin 
averaged (i.e., averaged over incident waves $1 and $2) and thermally averaged change in flux as 
a function of tip position, A/i(i?t, E) = fi{Rt, E) — /i(-Rt = 00, E), is given by: 

' m* \ 12 \Ep J J ' irup,dy 

2e h ^cos(9t) + cos(9f)- - ,. ^ , , , 

X sinh"^ [2XT{rtip,d - rk,d)] g{'2^{rtip,d - rk,d),T, Ep) + h.c. 




7r2 m" \ L2 \ lip / I ^ [rtip^d 
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_2 cos^^ + F{Rt) exp {i2kprup,d) + h.c. (56) 



where F{Rt) is some nonoscillatory function which depends on the particular configuration of 
scatterers, along with a some angular dependence of the tip and a exponential damping factor 
depending upon the tip position from the scatterers. Spin-orbit effects are not seen in Eq. (l56t due 
to the fact that for electrons moving along effective one-dimensional trajectories, the amount of 
spin rotation induced is simply proportional to the net distance the electrons have traversed. Since 
an electron which is directly scattered back to the detector has effectively traveled no net distance, 
the net spin rotation is zero. Another way to see this is that U{kark,d, 0'^)U{kare,k, ^e) — ^ when 
the emitter and detector are one and the same, since = r^.^^ and 6^^ = 6'g + vr. As shown at the 
end of Appendix C, the above conclusions also hold if the approximation to Gq{R) in Eq. ( l39t is 
not made. 

Consider the case explicitly illustrated in Figure El where the detector and the emitter are two 
distinct entities. In this case, an electron does not traverse the same path back to the detector, 
and the effects of spin-orbit coupling do not average away in the flux calculation. Performing the 
thermal averaging and spin averaging [using Eqs. (I45t-(l46bl. Afi{Rt) can be written as: 



m* I 12 \ Ef J I * re,tiprup,d 



't , 



+ exp(i/ci7'rs) exp(ifcc,r5) + G'2(^t) exp(-iA;ar5) 

+ ex^iikprs) {G^{Rt) exp(iA;„rD) + Gi{Rt) exp(-iA;„r£,) j + h.c. (57) 

where rg = r^^up + rtip,d and r^, = rup^d - re,up, and where Gi{Rt), G2{Rt), G^iRt), and G^iRt) 
are nonoscillatory functions of tip position which depend upon the particular configuration of 
scatterers. From Eq. (l57b . the expected elliptical fringes spaced at kpVs = 2nn with the detec- 
tor and the emitter acting as the foci of the ellipse are present; however, the amplitude of these 
oscillations are now modulated by the Rashba spin-orbit coupling. Since the electron trajectories 
from the emitter to the detector QPC are now two-dimensional, the electron's spin will undergo 
a trajectory-dependent spin rotation for each pathway between the emitter and the detector QPC. 
The interference between different pathways will thus have an additional, spin-dependent modu- 
lation. Such an amplitude modulation is similar to the Elliot- Yafet model of spin dephasing in the 
presence of spin-orbit coupling^*^. The first two oscillatory terms in Eq. dSTt . Gi{Rt) and G2{Rt), 
lead to elliptical amplitude modulations of the regular fringes spaced at kaf's = 2nn, while the G3 
and G4 terms in Eq. (BTt lead to a hyperbolic amplitude modulation spaced at kaTo = 2n7r, with 
the foci of the hyperbola again being the detector and the emitter. Interference between the terms 
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Gi and Gs and between the terms G2 and G4 lead to fringes at kare,tip = 2n7r for rup^d =constant. 
Likewise, interference between the terms Gi and G4, and between the terms G2 and G3 lead to 
fringes at kaVup^d = Snvr for r^^tip =constant. The presence of all four types of fringe patterns 
leads to a checkered pattern in A/i(/?f). This can be seen in the calculation of A/i(_R() in Eq. (1571) 
which is shown in Fig. |3l(A). In this simulation, the emitter was placed at = 1.5/i mX and the 
detector was placed Sfim away at = —l.bfimX . In addition, the thermal average of Afi{Rt) 
was evaluated by numerically integrating Eq. (l52t over the interval Ep ± GksT, and the following 
parameters were used (similar to the parameters found for InAs^^): m* = 0.022mo (where uiq is 
the free electron mass), Ep = 16 meV, T = 3K, and a = 4 x 10~^^eVm which gives a spin rota- 
tion length (i.e., the length required to rotate the spin by 180° of 1^, = nfi^ / (2m*a)) of 134 nm (For 
comparison, spin rotation lengths of ~ 1.8/um were found for heterostructures of GaAs/AlGaAs 
in a past study^^). Scatterers were randomly placed in the region [Y, X] =(0 fim, 6 /im)x (— 6 /^m, 
6 /im) with a scatterer density of 20 scatterers per /im^. All scatterers were modeled as cylindrical 
wells or barriers, with the well depths randomly chosen between ±.04 eV. The coupling constants, 
tk, were evaluated using Eq. (l37t and the results in Appendix A. The radius of the barriers/wells 
were all taken to be a = 3 nm, which gave kpa ^ .3 so that a model of "s"-wave scatterers 
could be used. The tip was modeled as a hard disc (i.e., infinite barrier) with the radius of the tip 
chosen to be a = 3 nm. Although the width of the actual depletion area induced by the tip in the 
2DEG is probably on the order of 100 nm, the above radius was chosen to be consistent with the 
"s"-wave model used in the calculations (further studies incorporating higher partial wave scatter- 
ing in the presence of Rashba coupling are currently being carried out and will be addressed at a 
later time). Figure|3tA) is a typical result from the calculations performed on numerous scatterer 
configurations. Besides the kinematical elliptical fringes spaced at kpVs = 2nn, the hyperbolic 
and elliptical modulations are clearly present in Figure |3] (A), along with the circular- like fringes 
about the emitter and the detector, which leads to a checkered pattern. 
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FIG. 3: Simulation of An{Rt) [Eq. d57t l for an unpolarized beam of electrons injected from an emitter 
located at re = 1.5 fimX and observed at a detector located at r ^ = —1.5 fj.mX (A) with and (B) without 
Rashba spin-orbit coupling for the same random configuration of scatterers. A scatterer density of 20 
scatterers per fim"^ was chosen; each scatterer was modeled as a cylindrical barrier/well of radius 3 nm with 
the height/depth of the potential randomly chosen between ±.04 eV. The tip was modeled as a hard disc 
of radius 3 nm. In (A), the amplitude of A/i(i?() is modulated by the Rashba spin-orbit interaction, due 
to the fact that electrons traveling from the emitter to the detector undergo net spin rotations. In (B), the 
expected elliptical fringes are observed. The same, arbitrary scale for Afi{Rt) was used in both (A) and 
(B). The following parameters were used in the simulation: m* = 0.022mo, T = 3 K, Ep = 16 meV, 
a = 4 X IQ-ii eVm. 
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For a comparison, simulations were also performed for the same scattering configurations and 
coupling constants but without the Rashba spin-orbit coupling (a = eVm). The fermi energy of 
these simulations, Ep, was chosen to be slightly higher in energy than Ep in Figure I^A) so that 
the magnitude of the fermi vectors was the same for both simulations, i.e., kp = kp. As expected, 
only regular elliptical fringes about the emitter and detector are shown in Figure |3tB), with any 
resulting modulation arising from the particular scatterer configuration. Note that the intensities 
of the fringes are larger when the tip is near to the detector than the corresponding fringes in the 
presence of Rashba coupling. Figure EtA) and (B) demonstrates this more clearly by plotting 
a slice of A/i(i?t) [shown in Figure |3tA)and (B)] along the F-axis, passing through the detector. 
Near the detector, the magnitude of Aij{Rt) is greater in the absence of Rashba coupling. However, 
the magnitudes of A/i(_Rt) with and without spin-orbit coupling are comparable far away from the 
detector. This is due to the fact that far away from the detector and emitter, the scattered wave 
functions can no longer resolve the detector and the emitter, i.e., U {kar^^d, d'^)U{kare,k, ) ~ 
which makes Afi{Rt) independent of a. It must be stressed, however, that although the form 
of the fringe pattern is robust to scatterer configurations, the overall intensity does depend on 
the scattering configuration. Figure EtC) and HD) give the same slice through a system with a 
different set of scatterers. 

B. Two Scatterer solution 

As mentioned earlier, if there is only one QPC, the effects of spin-orbit coupling are not ob- 
served if the electron trajectories from and towards the QPC are purely one-dimensional. As 
shown in Figure|5l two-dimensional multiple- scattering trajectories exist for electrons exiting and 
arriving at the detector. Consider the case of two scatterers: a moveable scatterer at Rt and a fixed 
scatterer at fg. In addition, both the tip and the fixed scatterer will be modeled as being infinite 
potential barriers of radius 3 nm. Using Eq. (l43b . the value of the wave function at each of the two 
scatterers is given by 



where A^^t = (1 —Det[Gl{rs)]) ^- This includes all orders of scattering between the two scatterers 
(i.e., any number of bounces between the two scatterers). The change in current as a function of 



*(r;) = Xs,Mrs) - Xs,tGl{fsMRt) 



(58) 



(59) 
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FIG. 4: Slices of Afi{Rt) through the detector at = — 1.5/im X. The plots shown in (A) and (B) come 
from Figure|3l Modulations are clearly evident in (A) due to spin-orbit coupling. Note also that the intensity 
of the fringes is larger in the absence of spin-orbit coupling [(B) vs. (A)]. However, the fringe intensities 
also depend upon the configuration of the random impurities. (C) and (D) show the same slice of Afi{Rt) 
for a different configuration of scatterers. Note again that modulations due to Rashba coupling are present 
in (C) and not in (D). 

Rt in the presence of the fixed scatterer at and a random configuration of weak scatterers can be 
found by inserting Eq. dSSt and Eq. (l59l) into Eq. dSOb and performing the thermal average using 
Eq. (l52b . Since t « 1 in the "s"-wave limit, it is useful to expand An{Rt) in powers of t. The 
single scattering contribution has already been discussed, and is given in Eq. (l56t . which is order 

—2 —3 

t . The next term of order t , which involves the interference between the trajectories shown in 
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FIG. 5: Possible multiple scattering trajectories which result in a net spin rotation applied to spins which 
make round trips from the emitter back to the emitter. Since the path is two-dimensional, a net spin rotation 
results, even when the detector and the emitter are the same. 

Fig. |5]and the single scattering trajectories, can be written as: 

A(3)/i(^t) = exp(iA;^s) {Ki{Rt) e^^^ik^Rs) + K2{Rt) e^V^^KRo)) 

+ exp(i^^5) (^7^3(^4) exp(-ifcc,.Rs) +i^'4(^t)exp(-iA;c.-RD)) 

+ exp(i^.R£,) {Li{Rt) ey^vi^KRs) + /^2(^t) exp(ifc„_R£,) j 

+ exp(iA;i?£)) {^^{Rt) exp(— I/cq^R^) + L^^Rt) exp{—ikaRD)j + h.c 

(60) 

where Rs = rup^d + Ts^tip and Rd = Uip^ - rs,tip- The functions Ki{Rt), K2{Rt), K^{,Rt), and 
Ki{Rt) (which mostly represent the interference between the impurity single scattering events 
and the trajectories shown in Fig. |5l) are nonoscillatory functions of Rt which depend upon the 
configuration of random scatterers, whereas the functions Li(Rt) , L2{Rt) , L^iRt) and L/^{Rt) 
(which represent the interference between the single scattering trajectories (for the tip and the 
fixed impurity) and the multiple-scattering trajectories shown in Fig. |5l) are nonoscillatory func- 
tions of Rt which only depend upon the position of the fixed scatterer at fg. Figure |6l shows a 
simulation of \^{^)^{Rt)\ in Eq. (1^ for a fixed, hard disc scatterer of radius 3 nm located at 
fs = 1.7 jjLmX + 2.4 jjivnY . The detector/emitter QPC was placed at (X, Y) = (0 /im, jj,m). Fig. 
^represents the evaluation of Eq. (l60b for the following parameters: T = 3 K, m* = 0.022mo, 

22 



a = 4 X 10~^^ eVm, and Ep = 16 meV (the same parameters as those given in Figure|3l). The ran- 
dom impurities were again modeled as cylindrical wells/barriers of radius 3 nm with depth/height 
randomly chosen between ±.04 eV and with a scatterer density of 20 scatterers per /xm^. Modula- 
tions due to spin-orbit coupling are again present, as predicted in Eq. (ISDt . 
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FIG. 6: Calculation of |A(3);u(^i)| in Eq. (I60t for a fixed hard disc scatterer of radius 3 nm located at 
fs = 1.7 fimX + 2.4 fimY. The modulations due to Rashba spin-orbit coupling are seen, similar to those 
shown in Figure |3l These modulations mainly result from the interference between the trajectories shown 
in Figure|51and the single scattering trajectories from the random impurities. The flux scale and parameters 



used in the calculation of |A(3)^(^t)| are the same as those given in Figure|3l 



If the scattering amplitudes, t, become large, then higher-orders (i.e., multiple bounces) must 
also be included. The interference between these different trajectories can lead to resonances 
induced by the scattering configuration. For the trajectories shown in Figure |5l however, no reso- 
nances due to spin rotation can be generated, since no net spin rotation is generated if the particle 
bounces from scatterer A to scatterer B and back to scatterer A again, as shown in Figure HA). 
The lack of spin rotation for such trajectories can be seen using the exact form of Gq{R) given in 
Eq. (l33t as follows: 



(61) 



However, for three or more scatterers(as shown in Figure EtB)), there exist trajectories which 
will give a net spin rotation. Calculating possible interference effects between multiple scattering 
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FIG. 7: Possible higher-order scattering processes. For two scatterers shown in (A), if the electron bounces 
between scatterer A and B, no net spin rotation results. For three or more scatterers, however, net spin 
rotation can occur. Such a possibility is shown in (B), where an electron traveling from A to C to B to A 
undergoes a spin rotation. 

trajectories requires using higher partial waves [Eq. (12^ 1 than the simple "s"-wave scattering 
models studied mostly in this paper, and will be investigated in the future. 

VI. CONCLUSIONS 

In this paper, a partial wave expansion for scattering from a cylindrical potential in the presence 
of Rashba spin-orbit coupling was developed and was used to construct an operator which gen- 
erates the scattered wave from the incident wave at scatterer k. This allowed for the development 
of point scattering models beyond the "s"-wave limit. The often studied "s"-wave scattering from 
delta-function potentials are shown to be different from the "s"-wave models developed in this 
work due to the fact that cylindrical wells/barriers do not scatter in the delta-function limit; how- 
ever, both models give the same qualitative results for the flux calculations presented in this work. 
Although only "s"-waves were discussed herein, extensions to higher partial wave scattering can 
be readily performed using the formalism developed in this work. Additionally, the operator 
and the calculated Green's function, G±{fi,f2, E), can be used to apply all the scattering theory 
machinery to study 2DEG confined in a variety of geometries. 

The Rashba spin-orbit coupling was shown to generate additional interference fringes in possi- 
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ble electron imaging experiments which were produced using a moveable scanning probe micro- 
scope tip. In the single-scattering limit, spin-orbit coupling doesn't produce any modulation in the 
observed flux using a single quantum point contact. This is due to the fact that no net spin rotation 
is generated from effective one-dimensional trajectories which start and end at the same location. 
If the injected current through a separate quantum point contact is measured instead, interference 
effects due to Rashba coupling are observed from the various two-dimensional trajectories from 
the emitter to the detector. This is due to the noncommutation of the resulting spin rotations along 
the trajectory, which results in spin-orbit-related interference effects. These interference effects 
are similar to the Elliott- Yafet mechanism of spin dephasing observed in electron systems. If 
multiple-scattering effects are also included, a single quantum point contact can again be used to 
observe the spin interference caused by the Rashba spin-orbit coupling. 

In the future, calculations involving higher partial waves and stronger scattering will be per- 
formed in order to look for possible spin resonances resulting from interference between the vari- 
ous multiple scattering trajectories. In addition, more realistic simulations of the scattering induced 
by a scanning probe microscope tip, requiring other partial waves in addition to the "s"-waves, will 
be performed. Finally, the multiple- scattering theory presented in this work can also be used to 
study scattering and polarization profiles generated in quantum wires where phase coherence ef- 
fects between the scattered waves can now be fully taken into account. 
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APPENDIX A: SOLUTION FOR SCATTERING FROM A BARRIERAVELL 



For the problem of the square well/barrier centered about ft, Eq. (QUI) and Eq. (Oil) require the 
various and df for |r — r/^ = a to satisfy the following equations: 

\fJi{k±a) + ^/hft^Hi{kia) + \/hfi^'^Hi{k2a) = ^/n^d^'^ Ji{Kia) + ./K^df^ Ji{K2a) 

±\fJi_i{k±a) + \/hfj^^Hi^i{kia) - >/k^fl^'^Hi_i{k2a) = ^/^Idf^ Ji^i{Kia) - y/i^df^ Ji^i{K2a) 

Xfk±j'i{k±a) + kl'^f^^H[{k^a) + kl'^f^^H[{k2a) = K^^^df^j[{K^a) + K^j^df^j[{K2a) 

±.k±\f Ji_^{k±a) + k^^"^ fj^^ Hi_^{kia) — k^^"^ f^'^Hi_^{k2a) = k^^"^ df^ Ji_^{Kia) — k^J"^ df^ Ji_i{K2a) 



where 

Af = 2i'exp(iA;± ■ ffc)exp(-i/6'o) (A2) 

Using Eq. (lAll) . the various values for the coefficients, df^, df^, f^^, and f^"^ can be found. In 
order to simplify the presentation of the solution to Eq. (lAll) . the following functions will be 
introduced to simplify the solutions: 

A/(p,g,a) = Ji{pa)Ji_i{qa) + Ji{qa)Ji^i{pa) 
AAi{p,q,a) = J'i{pa)J'i_i{qa) + J'i{qa)J'i_^{pa) 

Slip, 1, a) = Ji-i{qa)Hi{pa) + {-if Ji{qa)Hi^i{pa) 

Fl'{k±,q,r,a) = {Ji^i{ra)Ji{k±a)±{-lfJi-i{k±a)Ji{ra)) 

- gAAf(g,r,a) {jUra)Jl{k^a) ± {-lfjUk^a)j;{ra)) (A3) 
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The solution to the above equations can be written as 

= A 
= A 



±1 _ ^± F^{k±, f^i, K2, a)Gf{k2, K2, Ki, a) - K2, Ki, a)G\{k2, ki, K2, a) 

' y/h{Gf{ki, Ki, K2, a)Gf{k2, K2, Ki, a) - G]{ki, K2, ki, a)G]{k2, ki, K2, a)) 

Til 
,± Ji 



^ki 

fr = A' 



±2 _ ^± Ki, K2, a)Gl{ki, H2, Hi, a) - F^{k±, K2, ^1, a)Gf{ki, /ti, /t2, a) 

' %.{G\{k2, Ki, K,2, a)Gl{ki, K2, ki, a) - Gf{k2, K2, ki, a)Gf{ki, ki, K2, a)) 



X 



r±2 



'ko 



^±1 _ ^^Ji-i{K2a)Ji{k±a) ± Ji{K2a)Ji^i{kj^a) 



' yK7A(Ki,K2,a) 

Af Fl{k±, Ki, K2, a){gf{ki, K2, a)Gf{k2, k,2, i^i, a) - g^{k2, K2)G}{ki, K2, ki, a)) 



+ 



7±2 _ .±Ji-iiK,ia)Ji{k±a) T Ji{Kia)Ji-i{k±a) 

Ij — Ai . r 

,/K2A(Ki,K2,a) 



A(ki, k,2, a) y/i^{G^{ki, Ki, ^2, a')G^{k2, k,2, 1^1, a) - G]{ki, k,2, 1^1, a)G}{k2, ki, K2, a)) 
Af Fi^jk^, /t2, f^i, a)igfiki, K2, a)G}{k2, /ti, K2, a) - g^{k2, K2)Gf{ki, Ki, K2, a)) 
A(ki, K2, a) y^{Gf{ki, Ki, K2, a)Gf{k2, ^2, ^i, a) - Gj{ki, K2, ki, a)G]{k2, ki, K2, a)) 



^ ><t F?ik±, 1^1, 1^2, a){gl{ki, Ki, a)Gf{k2, K2, /ti, a) + g'^{k2, Ki)Gj{ki, K2, Ki, a)) 
A(ki, K2, a) ^/K^{Gf{ki, Ki, K2, a)Gf{k2, ^2, ^i, a) - G}{ki, K2, ki, a)G]{k2, ki, K2, a)) 

Xf Fl{k±, K2, Hi, a){gl{ki, ni, a)G]{k2, ft^i, /t2, a) + g'^{k2, Ki)Gf{ki, ni, K2, a)) 

A(ki, K2, a) ^/K^^Gf^ki, ki, K2, a)Gf{k2, K2, ki, a) - G\{ki, K2, ki, a)G\{k2, ki, K2, a)) 

(A4) 

It is useful to consider the limiting case of a hard disc, i.e., Vq ^ 00. In this case, 



m*a . /2m* I Vol 



m*a . 2m*\Vo\ ^ ^ 

K2 = — — +iW A (A5) 



with I Vol 00. In this limit, the coefficients are given by: 

(.±1 _ A;=^ Ji{k±a)Hi_i{k2a) ± Ji-i{k±a)Hi{k2a) 



fi- 



ki Hi{kia)Hi_i{k2a) + Hi{k2a)Hi_i{kia) 



r±2 Ji{k±a)Hi_i{kia)T Ji-i{k±a)Hi{kia) 

J I = 77=7771 n \ , tt n V77 n V (^6) 



^ Hi{kia)Hi^i{k2a) + Hi{k2a)Hi^i{kia) 

Note also that in the opposite limit of an infinite well, Vq ^ — 00, the expressions for fj^^ and fj^'^ 
are the same as those given in Eq 
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APPENDIX B: CONSTRUCTION OF Di 



For an arbitrary plane wave state specifed by energy E and with incident momentum vectors 
making an angle of 9q with respect to the F-axis, 

an operator Di can be constructed such that 



fl=l44--i^l (B3) 



First define the operator Pi by 



|/| dRx ORy 

with Pq = 1. Exponential functions of the form exp(ifc ■ R), where k = k{cos{9o)Y + sin(6'o)X), 
are eigenf unctions of Pi, where 

Piexp{ik-R) = k^^Uxp{iWo) exp{ik ■ R) (B4) 

The operator Di can be decomposed in terms of the operators Pi as follows: 

/ aiPi biPi+i\ 
A = i (B5) 
\ciPi^i diPi J 

where the coefficients, ai,bi, q and di need to be determined. Operating Di on \E'(i?) in Eq. (|Blll : 

'^\R)\ _ Aexp{iWo)exp{ih-R) f aikf^ + bik^+^^ \ 

Bexp{iWo) exp{ik2 ■ R) f o-ikf - 6; 4'^^' 

^2 \- exp{-i9Q){dikl^ - cik\ 



D 



^ ^ \- exp{-ieQ){dikl^ - cikl ^' 

, , f^HR)\ 

exp(i/^o) ^ (B6) 



The various coefficients must therefore satisfy 

aik^^ + fo/Zcf 
aikf - bik^2^^^ 
cikf'^^ + dikl^ 
ii 
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which gives 



ai 



hi 



Cl 





"-2 


+ kl^'\ 




^'\kl 






k\ 






^'\kl 






k\ 




k\ 


-'\kl 

"-2 


~|~ k-^ 
+ A:r^l 


^2 


-'^kl 


~|~ k-^ k<2 



For an arbitrary eigenstate of Hq [Eq. ©J with energy E, 



^^{R)) Jo \ V2 \exp{-ieo)J V2 V-exp(-ieo; 

dOo^ieo) (B9) 







Operating Di upon \l/(-R) gives 

^'Uk^V ^ I '^^o^^P(^^^o)'^^^o) ^^^^^ 

APPENDIX C: THE GREEN'S FUNCTION IN THE PRESENCE OF SPIN-ORBIT COUPLING 

The Green's function in the presence of Rashba spin-orbit interaction has a simple form in 
momentum space and can be written as 

G±(E) = -Km ^ 



'-^^H -E±ie 

1 ^* f ,r \^){k\ ^ , . ^ , . 

= - hm / dk— — — (l + cTx cos(0g) - cry sm(0^^ 

— lim - — — ■ / dA; = ( 1 — a v- cosfc/)^) + dy smichj:) ) 



(Cl) 



where (/)^ is the angle k makes with the respect to the F-axis. The form of the Green's function 
in Eq. (ICll) has been used in numerous studies of spin dynamics in 2DEGS. The position space 
representation of the Green's function, G±{ri,r2, E) = {ri\G±{E)\r2), however, has not been 
used to the best of the authors' knowledge. In the presence of a scattering potential, V{r), the total 
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wave function with energy E is related to G±{ri,f2,E) and V{r) by the Lipmann-Schwinger 
equation: 

^(^) = <l>(^) + j dfG+{R,f,E)V{r)<il{r) (C2) 

— * 

where $(-R) would be the wave function in the absence of the scattering potential V{r). 

Before calculating G±{fi, •r2,E), it is worth noting that the Rashba Hamiltonian, Hq [Eq. (jSl)], 
is invariant to combined rotations in spin and space about the z axis: 

^0 = M{e)HoM\e) (C3) 
M{e) = exp (^-i^Lz^ exp (^-i^^r^^ 

= R{e)z{e) (C4) 

Due to the above symmetry, the Green's function is also invariant: 

G±{E) = M{e)G±{E)M\e) (C5) 

From Eq. (IC5I) . it follows that 

{n\G±{E)\f2) = {n\M^e)M{e)G±{E)M^{e)M{e)\f2) 

= Z\9){R{9)n\G±{E)\R{9)r2)Z{9) (C6) 



G±{fi,f2, E) can be written as 



-m 



exp [k ■ r]l exp [k ■ f] 1 



G±{fi,f2-,E) = lim- — — -r / dk — „ , ,, h ^ ^ 

+ exD (k f] I ^QsC^fc) - Mh) _ cos(0g) - ay sin(0g) \ 

where r = ri — r2. In evaluating Eq. (IC7I) . it is advantageous to take r to be along say the F-axis. 
If f = r[cos(6')F + sin(6')X], then R{9)\f) = \rY). Therefore, using Eq. Ci), Eq. CH) can be 
written as: 

,^ ^ ^, — m* Z'^'" f°° I kexp(ikr cos((f)))l A; exp (iA;r cos (</))) 1 
G±(fi,r2,E) = hm— — - / d(P dk[ , . , + ^^^^ 



o{27Thy Jo ^ Jo \{k-ki±ie){k + k2) {k + ki){k - k2 ±ie) 

^ /Z^{9){axcos{(j)) -aYsm{(j)))Z{9) _ Z^{9){ax cos{(j)) - dy sm{(j)))Z{9)\ 
V {k-ki){k + k2±k) {k + ki±ie){k-k2) ; 

(C8) 



30 



where 



m*a I m*a\ 2m* E 



rrfa f m*a\ 2m/' E 



The integrals in Eq. dCSI) can be readily evaluated. The term in G±{fi,f2, E) proportional to 
the identity matrix is given by: 

— m* /"^'^ Ai ( k exp{ikr cos{(j))) k exp{ikr cos{(f))) \^ 

^(2^70 Jo \{k-ki±ie){k + k2) ^ {k + ki){k-k2±ie)) 

= lim^ TdA^^M^ ( I I I \l 

e^o 271%^ Jo ki + k2 \k + ki k - ki±ie k + k2 k - k2±ie J 

which is just the weighted average of two free particle Green's function with different wave vec- 
tors, ki and k2. 

Using the following two integrals: 

r-27r 

d</)exp(ifcr cos((/))) sin((/)) = 

d</)exp(iA;r cos(0))A; cos(0) = —Itt / d</)exp(i/cr cos((/))) 
or Jq 

d 

= -i^27rJo(fcr) 



(Cll) 



The terms in Eq. (IC8I) proportional to ax and ay are given by 



^o{2nh)^Jo Jo {k-ki){k + k2±ie) 

kZ\9){dx cos(0) - dy sin(0))Z(e) 



{k + k^±ie){k-k2) 
m*i - ^, a r Jo{kr) (I 1 1 1 



^^0 27ih^ dr Jo ki + k2\k + ki 

-Z\6)dxZ{6)m* d 



2{ki + k2)h^ dr 



ki + k2 \k + ki k — ki ± ie k + k2 /c — A;2 ± ie 
{H^ik^r) - H^{k2r)) 
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The Green's function, G±{ri,r2, E), can finally be written as 

m" ki ( H^{kir) iiff (fcir) exp (16 



G±{n,r2,E) = -1 



m 

1 



2^^ ki + k2 y iHf{kir) exp(-10) H^{kir) 

k2 I H^ihr) -iH^{k2r)exp{i6 



2r ki + k2 \ -iHf{k2r) exp(-10) H^{k2r) 



(C13) 



As noted earlier, G±(ri, r2, E) Is similar In form to Gq{R) operator [Eq. (TOl l found In the partial 
wave scattering analysis given In the text. It Is worth pointing out that, as was the case for Gq(R), 
G±{fi, r2, E)G±{f2, ri,E) oc 1, I.e., no net spin rotation Is observed when the particle traverses 
no net distance along a one-dlmenslonal path. 

Consider the experiment shown In Figure |2l for the case when the detector and the emitter are 
one and the same. With the approximation to Gq{R) made In Eq. ( l39b . It was concluded that 
Rashba spln-orblt coupling does not give rise to modulations In the net current as a function of Rf 
[Eq. (15^ 1. It Is worth pointing out that the conclusion reached by Eq. (15^ Is valid, even If the 
approximation In Eq. (l39t Is not made. Treating the emitter as a point source, the scattered wave 
function at the detector/emitter Is given by Eq. (l32t : 

N 



k=l 
N 



= J2G',iu)G+ifk,r,,E)ri (C14) 

k=l 

where rj represents the spin state of the Injected electron (can be taken to be either (p) or d") or 
some linear combination of the two). It Is easy to show, however, that GQ{fd)G+{fk, fa, E) Is not 
proportional to the Identity matrix. For kvk^ 3> 1, G'^{rd)G +{rk^ r^, E) can be written as 

e;(f,)e44,r.,£) oc a. 6.expK')^| ^^^^^ 

^Kd \6fcexp(-l^^) Ofc 



where 




^2,1 , [h 



kl ,2 / ^2 ,1 
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(C16) 



Eq. (IC15I) indicates that the incident spin state, rf, is modified after following the effective ID 
trajectory (although Eq. (IC15I) doesn't represent a spin rotation since it is not unitary). However, 
the only distance dependent factor in Eq. dClSI I contains krk,d and does not depend upon ka. 
Calculations of Afi{Rt) therefore will not possess any modulations which depend upon the Rashba 
interaction, ka, supporting the original conclusions made in Eq. (l56b . 
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